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Abstract. We present preliminary measurements of the bj and .fT-band 
luminosity functions (LFs) of field galaxies obtained from optical and K- 
band imaging of a sample of galaxies selected from the Stromlo-APM 
Redshift Survey. The bj LF is consistent with that previously published 
from photographic data. The K-band LF has been estimated over a range 
of 12 magnitudes and is reasonably well fit by a Schechter function with 
faint-end slope a = —1.2. 



1. Introduction 

Deep, near-infrared JC-band (2.2/im) galaxy surveys are a powerful tool for 
studying galaxy evolution (eg. Gardner et al. 1993, Cowie et al. 1994, Glaze- 
brook et al. 1995). Compared to blue-optical light, near-infrared light is a better 
tracer of mass in evolved stars and the correction for redshift dimming (the "fc- 
correction") is approximately independent of morphological type. The rapid 
evolution in galaxy luminosity apparent in the bj band is not seen in the K 
band. However, it is vital to have a reliable determination of the X-band lu- 
minosity function for nearby galaxies in order to interpret faint galaxy counts 
and to calculate the clustering of -fT-selected galaxy samples. The largest lo- 
cal -KT-band sample of galaxies with redshifts is that of Gardner et al. (1997). 
They measured redshifts for 510 galaxies selected from a -KT-band limited survey 
covering 4.4 square degrees. Since their survey is flux-limited, the majority of 
galaxies have if -band luminosities close to L* K . Thus they are able to measure 
the if -band luminosity function over a range of only 5 magnitudes, and the 
faint-end slope of their best-fit Schechter function, so important for predicting 
galaxy number counts, is poorly constrained (their Figure 1). 

One can improve on current estimates of the if -band LF without a huge 
investment of telescope time by observing galaxies selected by their intrinsic 
luminosity rather than their apparent flux. The Stromlo-APM galaxy survey 
(Loveday et al. 1996) is an ideal source for a new determination of the joint 
optical/near- infrared luminosity function ^(L^j, Lk) since redshifts have already 
been measured for 1797 galaxies with bj < 17.15 over a very large volume of 
space. The solid angle of the survey is 1.3 sr and the median redshift is about 
15,300 km/s. One can make use of the fact that K and bj luminosities are 
correlated, so that we can preferentially select galaxies of high and low luminosity 
and thus sample the luminosity range more evenly than a flux-limited sample. 



We are thus able to measure the luminosity function to fainter luminosities than 
from a flux-limited sample of similar size. 

2. Sample Selection 

Our aim in selecting a subset of Stromlo-APM galaxies for which ta obtain K- 
band photometry was to sample the magnitude range —22 < 13 < —13 as 
uniformly as possible. An added complication in defining the sample arose be- 
cause we wished to obtain optical CCD images for the same sample of galaxies. 
One planned use of this optical imaging is to measure morphological parameters 
for a representative sample of galaxies at low redshift in order to compare with 
HST observations of galaxies at high redshift (z > 0.4). To obtain compara- 
ble linear resolution to the HST data required observing galaxies at z < 0.04, 
assuming ground-based seeing of 1.3 arcsecond. Our "primary" sample thus 
consists of galaxies at redshifts z < 0.04. We divided the magnitude interval 
—22 < Mbj < — 13 into 90 bins each of width 0.1 mag. We then randomly 
selected up to six galaxies from the Stromlo-APM survey with z < 0.04 in each 
bin. Due to its redshift limit of z < 0.04, this primary sample contains rather 
few galaxies brighter than Mbj = —20. We therefore formed a supplementary 
sample, consisting of galaxies at z > 0.04 to "top up" each magnitude bin, 
where possible, to six galaxies. This supplementary sample consists entirely of 
galaxies with Mf,j < —20. The primary sample contains 283 galaxies, and the 
supplementary sample contains 80 galaxies, giving a total sample size of 363 
galaxies. 

3. Observations 

.ftT-band imaging of the above sample of galaxies was carried out at the Cerro 
Tololo Interamerican Observatory (CTIO) 1.5m telescope using the CIRIM in- 
frared array over the nine nights 1996 August 31 - September 4 and 1997 October 
19-22. The pixel size at f/7.5 is 1.16", allowing most galaxies to be observed at 
9 non-overlapping positions on the 256 x 256 array. Total integration time for 
each galaxy was 300s. The infrafred frames were reduced using IRAF, and pho- 
tometry was performed using SExtractor 2.0.8 (Bertin & Arnouts 1996) with the 
"mag_best" option. Magnitude errors were estimated by combining in quadra- 
ture SExtractor's estimate of the error from photon statistics and the difference 
between magnitudes measured using local and global estimates of the sky back- 
ground. 343 galaxies were observed with an estimated K-band magnitude error 
of less than 0.3 mag (rms mag error = 0.06 mag) and were calibrated using 
standard stars observed from the list of Elias et al. (1982). 

Figure [l] shows the rest- frame (bj — k) versus Mk colour-magnitude relation 
for our data. The fit to all galaxies is given by 

(bj - k) = -0.260 x M K - 2.11, a = 0.86. (1) 

We assume a -KT-band k-correction of — 2.5z for all galaxy types. 



1 Throughout, we assume a Hubble constant of Ho = 100 km/s/Mpc. 
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Figure 1. bj — k versus colour- magnitude plot using APM bj 
magnitudes. Plus signs represent early-type galaxies, asterisks late- 
type galaxies and dots represent unclassified galaxies. The solid line 
shows a least-squares fit to all galaxies, the dahsed line a fit to early 
types and the dot-dashed line a fit to late type galaxies. 



Optical imaging was performed in the U, B and R bands using the CTIO 
1.5m telescope with a Tex 2048 2 CCD over the ten nights 1996 September 7- 
16. Integration times were 120s in R, 240s in B and 120s in U. The [/-band 
exposures were too short to provide accurate galaxy photometry but were taken 
under photometric conditions, allowing the possibility of later calibration of 
deeper, non-photometric [/-band observations. Galaxy photometry was done 
using an earlier beta-release (1.2b9b) of SExtractor and so the photometry pre- 
sented here is only preliminary. 300 galaxies have reliable B and R magnitudes, 
and are calibrated with Landolt (1992) standards. The colour equations of Couch 
& Newell (1980) were used to obtain a bj magnitude from B and R. Figure |2| 
plots these CCD bj magnitudes against APM bj magnitudes. The mean and 
rms APM — CCD magnitude is Am = 0.09 ± 1.05. This scatter is larger than 
the 0.3 mag scatter for the full, m > 15 Stromlo-APM sample (Loveday et al. 
1992) since 1) galaxies with APM magnitude brighter than 15 are included in 
the current sample, these are badly saturated on the photographic plates, and 
2) the preferential sampling of galaxies of very high and low luminosity means 
that galaxies with poor APM magnitudes are more likely to be included. For 
example, the outlying galaxies in the lower-right of Figure |2| are due to the APM 
machine measuring just a small part of a large spiral galaxy, and thus grossly 
underestimating the galaxy's luminosity. 
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4. Luminosity Functions 

When one has a sample selected on one quantity mi (in our case APM bj 
magnitude) and wishes to estimate the luminosity function for another quantity 
777,2 (eg- CCD bj or K magnitude), the best way to proceed is to calculate a 
bivariate luminosity function (BLF) <j)(Li,L-2) allowing for any selection effects 
in L\ and then to integrate over L\ to obtain 4>(Li ). One can estimate the shape 
of (p(Li,L2), independently of inhomogeneities in the galaxy distribution using 
the maximum likelihood method of Sandage, Tamman and Yahil (1979). The 
probability of seeing a galaxy with luminosities L\ and L\ at redshift zi is given 
by 

max 

/ / (f)(Li, L2)S(Li)dLidL 2 . (2) 

^2 min m i n (Zi) 

The function S(L\) accounts for the known selection in L\ and the luminosity 
limits -kimin^i) an d iimax(^) are the minimum and maximum luminosities 
observable at redshift Zi in a sample limited by apparent m\ magnitude. If there 
are no flux limits in the 7772-band, then the integral over L 2 runs from to +oo. 
The maximum-likelihood shape of the BLF 4>{L\,L2) is estimated by maximizing 

the likelihood C = Ylfji Pi (the product of the individual probabilities pi for the 
N g galaxies in the sample) with respect to the parameters describing the BLF. 

In practice, we do not have a good a priori parametric model for (f>(Li, L 2 ), 
and so instead we measure 4>(Li,Li2) in a non-parametric way using an exten- 
sion of the Efstathiou, Ellis and Peterson (1988) stepwise maximum likelihood 
(SWML) method. Sodre and Lahav (1993) have extended the SWML method 
to estimate the bivariate diameter-luminosity function and to allow for sam- 
ple incompleteness. We adopt their extension of the SWML estimator here, 
including the sampling function S(L\) separately for the primary and supple- 
mentary galaxy samples. We normalise our LFs to the mean density of galaxies 
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Figure 3. The iT-band luminosity function estimated from nine 
Soneira- Peebles simulations (points) along with the input LF (curve). 



with -22 < M bj < -13 in the full Stromlo-APM sample, n = 0.071/i 3 Mpc~ 3 , 
calculated as described by Loveday et al. (1992). 

Once one has obtained the SWML estimate of </>(-Li, L2), one can integrate 
over L\ to obtain <j)(L2 ) and then fit a given functional form, eg. a Schechter 
(1976) function, by least-squares. 

4.1. Test of the Method 

We have tested the above procedure by using it to estimate the -fT-band lu- 
minosity function from a set of Monte Carlo simulations. We generated nine 
mock Stromlo surveys by a Soneira and Peebles (1978) hierarchical clustering 
simulation. Each galaxy in the simulation was assigned a K-band luminosity 
drawn at random from a Schechter function with a = —1.21 and = —24.7. 
Each galaxy was then assigned a bj magnitude according to our observed colour- 
luminosity relation (Q). Galaxies were selected on their apparent bj magnitude, 
bj < 17.15. This process was repeated until each simulation contained 2000 
galaxies. We then sampled each simulation by absolute M bj magnitude as de- 
scribed in §^], finally yielding an average of 359 galaxies per simulation. We 
calculated the K-band luminosity function (f)(Lx) for each simulation as de- 
scribed in and fit a Schechter function to each by least squares. Averaging 
over the nine simulations, and estimating the BLF in bins of width 0.5 mag, 
we measure mean and rms Schechter function parameters a = —1.16 ± 0.06, 
M* = —24.4 ±0.6. The errors on the mean values are \/9 times smaller than the 
quoted rms scatter between the simulations. Thus our estimate of M* is biased 
1.5cr too faint and a is overestimated (too shallow) by about 2.5o\ However, our 
estimates are within the la error from a single realisation. 
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Figure 4. The bj luminosity function estimated from our CCD mag- 
nitudes (filled symbols). Also shown (dashed line) is the Schechter 
function fit by Loveday et al. (1992) from photographic magnitudes. 



The SWML estimates of the K-band LF for the simulations are shown 
in Figure |3[ The points show the mean SWML estimate of 4>(Lk) from the 
nine simulations and the curve shows the input Schechter function with shape 
a = —1.21, = —24.7. The error bars going through the data points show 
the rms scatter between realisations. The error bars offset slightly to the right 
show the mean error predicted by the covariance matrix (see Efstathiou et al) 
and are in reasonable agreement with the rms scatter between realisations. 

Overall, we find that our procedure for estimating (\){Lk) from a sample 
limited by apparent b magnitude and futher selected by absolute B magnitude 
provides a robust and only weakly biased estimate of the /T-band LF over a wide 
range of absolute magnitudes. 

4.2. <P(L bj ) 

Our estimated 6j-band luminosity function is shown in Figure [|. The curve 
shows a Schechter function fit to the SWML estimate using least squares and 
allowing for finite bin width. The best-fit Schechter parameters are a = —1.16, 
Ml = -19.52 and (f>* = 0.018/i 3 Mpc~ 3 . Also shown in this figure is the 
Schechter function fit by Loveday et al. (1992) from photographic magnitudes. 
Despite the large scatter in our CCD versus APM bj magnitudes (Fig. |2|), we 
find that the two estimates of the shape of 4>(Mbj) are in reasonable agreement. 
Although the new estimate has slightly steeper faint-end slope, the slopes are in 
fact consistent within the la uncertainties. 
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Figure 5. The K-band luminosity function estimated from our sam- 
ple (filled symbols). Also shown are the results of Gardner et al. (open 
circles) and Szokoly et al. (open squares). 



4.3. <f>(L K ) 

Our estimated .fT-band luminosity function is shown in Figure [j| The curve 
shows a Schechter function fit to the SWML estimate using least squares and 
allowing for finite bin width. The best-fit Schechter parameters are a = —1.22, 
Mk = —24.73 and 0* = 0.0073/i 3 Mpc~ 3 . Also shown in this figure are recent 
estimates of the /f-band LF from the X-selected samples of Gardner et al. (1997) 
and of Szokoly et al. (1998). These estimates are in good agreement with ours 
at Mk ~ —22 but fall below our new estimate at the bright end. The reason 
for this is almost certainly the fact that our sample is selected in bj, whereas 
the Gardner et al. and Szokoly et al. samples are -fT-selected. In a blue-selected 
sample, the bright end of the K-band LF will be dominated by red galaxies and 
the faint end dominated by blue galaxies. Since red galaxies tend to be luminous 
ellipticals, this would explain why we see a larger 4>{Mk) at the bright end than 
seen in K-selected samples. 

5. Conclusions 

We have presented preliminary estimates of the bj and K-band luminosity func- 
tions obtained from CCD and infrared array imaging. Our bj LF is consistent 
with an earlier estimate from photographic data. We measure 4>(Lk) over a 
range of 12 magnitudes, a significantly greater range of luminosities than has 
been measured until now. We are thus able to place much tighter constraints on 
the faint-end slope of the if-band LF, albeit for a 6j-selected sample. Planned 
future work includes measurement of the bivariate LF cf)(Lb r , Lk) using CCD bj 
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magnitudes, the .R-band LF and luminosity functions of galaxies subdivided by 
colour and morphological type. Given the topic of this meeting, it will also be of 
great interest to calculate bivariate luminosity-surface brightness distributions. 
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